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Abstract 

Harboring a large number of endemic species, tlie Alps and the Western Carpathians are considered as major centers of 
biodiversity. Nonetheless, the general opinion until the turn of the millennium was that both Central European mountain 
regions did not provide suitable habitat during the Last Glacial IVlaximum, but were colonized later from southern refuges. 
However, recent molecular genetic studies provide new evidence for peripheral Alpine refuges. We studied the 
phylogeography of the calciphilous land snail O. dolium across its distribution in the Alps and the Western Carpathians to 
assess the amount of intraspecific differentiation and to detect potential glacial refuges. A partial sequence of the 
mitochondrial COI was analyzed in 373 specimens from 135 sampling sites, and for a subset of individuals, partial sequences 
of the mitochondrial 16S and the nuclear histone H3 and H4 were sequenced. A molecular clock analysis was combined 
with a reconstruction of the species' geographic range history to estimate how its lineages spread in the course of time. In 
order to obtain further information on the species' past distribution, we also screened its extensive Pleistocene fossil record. 
The reconstruction of geographic range history suggests that O. dolium is of Western Carpathian origin and diversified 
already around the Miocene-Pliocene boundary. The fossil record supports the species' presence at more than 40 sites 
during the last glacial and earlier cold periods, most of them in the Western Carpathians and the Pannonian Basin. The 
populations of O. dolium display a high genetic diversity with maximum intraspecific p-distances of 18.4% (COI) and 14.4% 
(16S).The existence of various diverged clades suggests the survival in several geographically separated refuges. Moreover, 
the sequence patterns provide evidence of multiple migrations between the Alps and the Western Carpathians. The results 
indicate that the Southern Calcareous Alps were probably colonized only during the Holocene. 
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Introduction 

The Pleistocene climate changes shaped the phylogeographic 
patterns of various organisms [1]. In particular, the severe cooling 
starting with the end of the Early Pleistocene (about 900 kya) was 
the starting point for massive glaciations in the northern 
hemisphere [2]. Mountainous regions such as the Central 
European Alps were heavily affected due to shifts in temperature 
and humidity, and the expansion of glaciers, resulting for many 
taxa in the fragmentation of populations, complete or local 
extinction, and the loss of variation due to genetic bottlenecks [1]. 
Thereby, the Last Glacial Maximum (LGM; 30-18 kya [3]) is of 
most relevance in respect to the current distribution of Central 
European species. Although the existence of glacial refuges at the 
periphery of the Alps was discussed more than half a century ago 
[4], the general opinion until end of the 20th century was that 
glacial refuges were located mainly in southern regions [1]. 
However, molecular genetic analyses and fossil data revealed the 
existence of northern refuges in the Western Carpathians [5,6] 
and the Pannonian Basin [7]. Several peripheral Alpine refuges 
were proposed for sihcophUous mountain plants [8] and 



calciphilous land snails such as Arianta arhustorum [9—12], Carychium 
minimum, Carychium tridentatum [13], Trochulus oreinos [14] and 
Trochulus villosus [15]. 

Most molecular genetic studies investigating glacial refuges were 
based on the assumption that populations diverged during 
isolation in geographically separated areas, and that populations 
of former refuge areas are now characterized by high genetic 
diversity and the presence of rare (private) alleles [16]. Neverthe- 
less, a major obstacle in identifying signatures of Pleistocene 
refuges is that phylogenetic signals are blurred because of 
migration and intermixture of previously separated populations. 
Therefore, species showing low dispersal capabilities and specific 
habitat requirements, which is the case in most land snails, might 
be suited to infer past distributional patterns. 

In the present study, we investigate the phylogeographic 
patterns of the land snaU Orcula dolium (Draparnaud, 1801). The 
species inhabits all major limestone areas of the Alps and the 
neighboring Western Carpathians [17,18]. 0. dolium is usually 
associated with mountainous forest habitats or rocky landscapes 
with patches of vegetation, but it is also found on rocky slopes at 
high altitudes up to 2160 m above sea level (asl) [17] (and data of 
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present study). Within the Alps, the east-west orientated Central 
Alps, consisting mainly of silicate rock, represent the largest 
distributional barrier for the species; they separate the Northern 
Calcareous Alps from the Southern Calcareous Alps. The Vienna 
Basin, closing the Pannonian Basin to the north, constitutes 
another distributional gap because; it s(;parates the populations of 
the Northern Calcareous Alps (including the Wienerwald) and the 
Western Carpathians. Recent investigations of loess profiles from 
the Pannonian Basin show that during the Late Pleistocene 0. 
dotium also occurred in the lowlands of the region [19], although 
the species seems to have vanished from it during the Holocene. 
The margins of the Northern Calcareous Alps and the Western 
Carpathians, enclosing the Vienna Basin east and west, harbor the 
majority of the described 23 subspecies [20]. As these areas were 
only partially glaciated during the LGM, they also come into 
consideration as potential rcfugc-s for the- species. 

We perform a comprehensive phylogcographic study of 0. 
dolium analyzing mt and nc markers to detect potential glacial 
refuges and to assess the amount of intraspecific variability. 
Furthermore, we test whether the described subspecies are 
dilferentiated genetically in the mitochondrial (mt) and nuclear 
(nc) markers. Our sample covers almost the entire range of the 
species. We include data of the Pleistocene fossil record to obtain 
insights into the species' past distribution. In order to estimate 
which potential areas were inhabited by ancestral populations of 
0. dolium and to trace the distribution patterns of the mt lineages 
throughout time, a molecular clock analysis is performed and 
combined with a phylogcographic range reconstruction. 

Methods 

Study Area and Sampling 

Specimens were collected from a large part of the species' 
distribution, including several Alpine and Western Carpathian 
type localities. 0. dolium is not protected by conservation laws of 
the countries where the collections were performed. Thus, in 
general, permissions were not necessary. For protected areas in 
Austria, permissions were provided by federal states authorities. 
Permit numbers: RU5-BE-64/01 1-2013 (Lower Austria), FA13C- 
53 Sch 6/6-2007 (Styria) and NlO-1 17-2008 (Upper Austria). 
Most samples of 0. dolium were collected in the Northern 
Calcareous Alps (Austria and Germany), a lesser fraction in the 
Western Carpathians (Slovakia), the Southern Calcareous Alps 
(Austria, Slovenia and Italy) and the Western Alps (Switzerland), 
totaling 373 specimens of 135 sites (Table 1). The habitats include 
wooded areas in the lowland, mountainous vegetated areas and 
rocky slopes in the alpine zones, with an altitudinal range from 
120 m to 2160 m asl. Elevation and position were determined via 
GPS. At every sample locality, if available, a minimum of three 
living specimens was collected, prepared for DNA analyses and 
stored in ?>()% ethanol following the protocol of [21]. Orcula conica 
(Rossmassler, 1837) (ID: 3899; Trogener Klamm, Carinthia) was 
used as outgroup. Selected type specimens of 0. dolium subspecies 
included in the present study are shown in Fig. SI. For the 
inference of substitution rates used in the molecular clock analysis, 
six specimens oiOrcuklla bulgarica (Hesse, 1915) and four specimens 
of Ouulella aragonica (VVesterlund, 1897) were included. Samples of 
Onuklla hulgarica were collected by Barna Pall-Gergely (Shinshu 
University, Matsumoto, Japan; Turkish samples) and Alexander 
Reischuetz (Greek samples). DNA samples of 0. aragonica were 
obtained from Benjamin Gomez-Moliner (Universidad del Pais 
Vasco, Vitoria, Spain). Voucher specimens of the first three taxa 
are deposited in the Natural History Museum Vienna, the 
whereabouts of the 0. aragonica vouchers are provided in [22]. 



Every individual sample was assigned an ID consisting of a unique 
specimen number and a locaUty tag. The latter encodes the Alpine 
geographic mountain region as classified in the SOIUSA system 
[23], with localities of each region numbered from west to east 
(Table 1). Due to their geographic vicinity to the Northern 
Calcareous Alps, the sites located in the Fischbach Alps, Eastern 
Styrian Prealps and Lavanttal Alps as well as those of the 
Wienerwald, which are classified to the Central Alps in the 
SOIUSA system, are treated as Northern Calcareous Alpine 
mountain regions here. A map providing an overview of the 
mountain regions in\'cstigatc'd is shown in Fig. 1. The Slovakian 
sites were eac:h assigned to one of the geological areas defined for 
the Carpathians [24]. To illustrate distribution patterns, the 
haplotypes in the phylogenetic trees and the histone network are 
marked by different colors, corresponding to the SOIUSA 
mountain regions as shown in Fig. 1. 

Molecular Markers and Primer Design 

DNA was extracted from 373 0. dolium specimens with the 
QIAgen Blood and Tissue Kit, using a piece of foot tissue. A 
partial region of the mt cytochrome c oxidase I gene (COI) was 
sequenced in all 373 specimens. From a subset of 54 individuals 
(including representatives of the major mt clades, type localities, 
peripheral geographic regions and contact zones of distinct mt 
clades), two additional markers were amplified: a section of the mt 
16S gene, as well as a nc sequence comprising almost the entire 
sequences of the histone 4 and 3 genes and the complete internal 
spacer region (H4/H3). The two histone genes are orientated in 
opposite direction and are separated by a non-coding spacer 
region, an arrangement which is probably not universal in 
gastropods but was verified so far particularly for species of the 
informal group of Orthurethra (sensu [25]). The COI and 16S 
sequences were also amplified in the eleven outgroup specimens. 
The COI forward primer COIfolmerFwd [14] is a \ ariant of the 
standard primer LCO1490 [26]; H2198-Alb [10] was used as 
reverse primer. New 16S primers were designed for the 
amplification of a fragment of approximately 850 bp. The forward 
primers, 16SL0rcl Jiu and 16SL0rc2Jiu, bind about 50 bp away 
from the 5'-end of the 16S rRNA gene and the reverse primer 
16SL0rc_rD is situated in a conserved region about 850 bp 
downstream. The H4/H3 primers, OrcH4_leftl and OrcH4_left2 
(positioned at the 3 '-end of the H4 gene) were designed based on 
alignments of H3 and H4 sequences from GenBank. The reverse 
primer OrcH3_rightl (at the 3'-end of the H3 gene) was published 
by [27]. Internal primers for sequencing [OrcH3S_right3 and 
OrcH4S_left3) were positioned close to the spacer in the coding 
sequences to obtain the complete 1100 bp fragment with t^vo 
sequencing runs. The PGR primers {OrcH4_kJil, OrcH4_left2 and 
OrcH3_rightl) cover a wider spectrum of Orthurethra taxa, while 
the internal primers {OrcH3S_right3 and OrcH4S_left3) were 
especieiUy adapted to the genus Orcula. All primers are listed in 
Table 2. 

PCR Amplification and Cloning 

COI and 16S fragments were amplified with the Roche Taq 
DNA polymerase for direct setjuencing. The PCR started with a 
denaturation step for 3 min at 94°C, followed by 35 (Ax les with 
30 s at 94°C, 30 s at the particular annealing temperature 
(Table 2), and 1 min at 72°C, followed by a final extension for 
7 min at 72°C. The PCR for the nc H4/H3 fragments was 
performed with the standard protocol of the Finnzymes Phusion 
polymerase, which has proofreading activity. PCR started with a 
denaturation step for 30 s at 98°C, followed by 35 cycles with 10 s 
at 98°C, 10 s at the particular annealing temperature (see Table 2), 
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Figure 1 . Distribution of mountain areas investigated in the present study. The colors correspond to those used In the mt trees (Figs. 1 and 
4) and in the histone H4/H3 network (Fig. 5). The outlines of the Northern Calcareous Alps (NCA), the Southern Calcareous Alps (SCA), the Western 
Alps (WA) and the Central Alps (CA) are framed in black. The names of the mountain regions and abbreviations are provided in the figure. 
doi:10.1371/journal.pone.0096012.g001 



and 30 s at 72°C, followed by a final extension for 7 min at 72°C. 
Purification and sequencing of the PGR products (both directions) 
were performed at LGC Genomics (Berlin, Germany), using the 
PGR primers for sequencing, except for the H4/H3 section, which 
was sequenced with the internal primers OrcH3S_right3 and 
OrcH4S_kft3 . The cloning of PGR products was performed for 
the H4/H3 primer design phase and for those individuals that 
proved to be heterozygous with respect to insertions and deletions 
(indels), resulting in varying spacer lengths and thus impeding 
direct sequencing. Nine of the 14 specimens with heterozygous 
sequences yielded fragments differing in spacer length. The 
fragments were purified using the QIAquick Gel Extraction Kit 
(QIAGEN), extended by A-endings with the DyNAzyme II DNA 
polymerase (Finnzymes) and cloned with the TOPO TA cloning 
kit (Invitrogen). Two to four clones were sequenced until the two 
main variants were obtained. In two cases (samples WIW2_392 
and MLF1_3939) more than two length variants were obtained, 
which is not completely unexpected in multi-copy gene families. 
Sequencing was also performed at EGG Genomics (Berlin) using 
Ml 3 universal primers. AH sequences are deposited in GenBank 
under the following accession numbers: KC568830-KG569204, 
KJ656162-KJ656172 (GOI); KG569205-KG569260, KJ6,56173- 
KJ656183 (16S); KC569261-KC569327 (H4/H3). 

Sequence Analysis and Phylogenetic Tree Calculation 

Sequences were edited using Bioedit 7.1.3 [28]. When directly 
sequenced H4/H3 fragments provided ambiguous positions, the 
respective sites were filled with the corresponding lUPAG codes. 
The complete COI data set comprises 374 sequences (including a 
single specimen of O. conica as outgroup). The alignment of the 655 
nucleotide sites was straightforward because there were no indels. 
Statistical analyses were performed using all GOI sequences. 
Haplotype and nucleotide diversity were calculated with DnaSP 



5.10 [29], and uncorrected genetic j6-distances between the clades 
of O. dolium and 0. conica were calculated with MEGA 5.2 [30] . For 
phylogenetic tree calculation, identical haplotypes from the same 
geographic areas (SOIUSA codes) were collapsed resulting in a 
total of 220 COI sequences of 0. dolium (197 haplotypes). Prior to 
the phylogenetic tree inference, a search for the best fitting 
substitution model was performed with jModelTest 0.1 [31]. A 
Bayesian inference (BI) was calculated with MrBayes 3.2.2 [32,33] 
for 5xl0'' generations (samplefreq= 100; nruns = 2; nchains = 4), 
applying the parameters obtained from the model test (GTR+G+I; 
nst = 6, rates = invgamma). Tracer 1.5 [34] was used to assess 
whether the two runs had converged and when the stationary 
phase was reached. The first 25% of the trees were discarded as 
burnin and a 50% majority rule consensus tree was calculated 
from the remaining trees. 

The 16S sequences (55 specimens including 0. conica) were 
aligned with GlustalX [35] using default parameters. The original 
alignment contained 879 positions of which all 61 gap positions 
were excluded with TrimAl 1.3 [36], implemented in the 
Phylemon 2.0 web tools [37], using the "no gap" option. Another 
84 sites were excluded by performing the "strict" option, leaving 
734 positions in the final alignment for the phylogenetic tree 
analyses. Of those, 217 sites were variable, compared to 331 in the 
original and 298 sites in the "gaps excluded" alignment. The "no 
gap" alignment was also used for calculation of uncorrected p- 
distances. The BI was performed with two data partitions (COI 
and 16S "strict"), using the substitution models suggested by 
jModelTest 0.1 [31] (COI: HKY+I+G: nst = 2, rates = invgamma; 
16S "strict": HKY+I+G: nst = 2, rates = invgamma), and allowing 
MrBayes to evaluate the model priors independentiy. A Maximum 
Likelihood (ML) tree was calculated with MEGA 5.2 [30], 
applying the sequence evolution model GTR+G (5 rate catego- 
ries)+I and performing 1000 bootstrap replicates with SPR 



PLOS ONE I www.plosone.org 



11 



May 2014 | Volume 9 | Issue 5 | e96012 



Phylogeography of the Land Snail Orcula dolium 



Table 2. Primer sets for amplification and sequencing. 




Region 


Primer (5' to 3') 


Origin 


Fragment size 


T Phusion 


T Roche 




CO! fwd 


CO!fotmer_fwd: 

/^y-~-T-/" A A A A T/" ATA A A AT A J 1 /-"/^ 

GG / C/l/\ LfKR TCA TAAA GA TA 7TG G 


[26] 


655 bp 


54X 


50°C 




CO! rev 


H2198-Alb: 

TA TA CTTCAGGA TGA CCAAAAAA TC 


[10] 










16S fwd 


leSiOrclJwd: 

1 1 A /•^TTTT/^ AT A AT^/^TTA A A f^T A 

TTALC TTTTG LA TAA / G G TTAAAC TA 


present 
study 


c. 850 bp 


60"C 


54°C 






16SLOrc2_fwd: 

I I A /"/ — J J J 1 A T A A Tf^ / — TT A A A TT A 

TTALLI 1 1 iGCA lAAlGGTTAAATTA 


present 
study 










16S rev 


16SL0rc_rev: 

/^/^^T/^T^ A A /^T/^ A ^ A T/" A T^ 


present 
study 










H4/H3 fwd 


OrcH4_leftl: 

GTGCGTCCCTGGCGCTTCA 


present 
study 


C.1100 bp 


7rc 


57°C 






OrcH4_left2: 

GGCCCTTCACCGCGTACACT 


present 
study 










H4/H3 rev 


OrcH3_nghtl: 

TGGGCATGATGGTGACACGCT 


[27] 










intern fwd 


OrcH4SJeft3: 

CGGTCTGAACTCAGATCATG 


present 
study 










intern rev 


OrcH3S_right3: 
CGGTCTGAACTCAGATCATG 


present 
study 











The annealing temperatures are provided for the Finnzymes Phusion and the RocheTaq Polymerase, respectively. 
doi:1 0.1 371 /journal.pone.009601 2.t002 



(Subtree-Pruning-Redrafting) as heuristic method for tree infer- 
ence. Support values of the ML analysis are provided for the 
combined COI and IBS tree. 

A Median-Joining network was calculated from the H4/H3 
data with Network 4.6.0.0 (Fluxus Technology Ltd.) using the 
default settings. Two networks were calculated: one without gaps 
and one keeping the gap positions. The networks were then post- 
processed to reduce unnecessary median vectors using the MP 
option. As the gap positions contained valuable information, we 
show the network including the gap sites. 

Molecular Clock Analysis and Geographic Range 
Reconstruction 

A reconstruction of the historical biogeography of 0. dolium was 
performed with Lagrange 2.0.1 [38]. Lagrange 2.0.1 uses a 
dispersal-extinction-cladogenesis (DEC) modeling, which allows 
analyzing the ML values of rate transitions as a function of time. 
The calculations were based on a molecular clock dated linearized 
BI tree calculated with BEAST 1.7.5 [39]. The fossil record of O. 
dolium comprises only material of the Holocene and the Middle 
and Late Pleistocene. This might be due to the fact that during 
these time periods climate conditions promoted the accumulation 
of loess sediments containing high numbers of gastropod shells. In 
contrast, fossil records of the Early Pleistocene and the Pliocene 
are very scarce and lack many land snail species endemic to 
Central Europe, among them 0. dolium. Consequently, dating the 
stem of the tree with the species first occurrence in the fossil record 
was not a reasonable option. Hence, substitution rates were 
inferred from a comparison of two other Orculidae, Orculella 
aragonita and 0. bulgarica. The latter is widespread in South-Eastern 
Europe and Western Asia, whereas its sister species, 0. aragonica, is 
distributed only in the Iberian Peninsula [22]. The earliest fossil 
record of 0. aragonica/ 0. bulgarica dates from about 1.8 mya at the 
Almenara-Casablanca karst complex (Castellon, Eastern Spain), 
which contains Miocene to Early Pleistocene sediments [40,41]. 
We assumed that the earliest record coincides with the time period 



when the ancestral lineages of the two species split from each 
other. Calculation of substitution rates and molecular clock 
analysis were both performed with the same COI and 16S 
("strict", 696 positions) alignments (including the sequences of 0. 
conica and the two Orculella species). For each partition, TN93+G 
was applied both to the measurement of substitution rates and the 
inference of the molecular clock dated tree. Molecular clock tests 
were performed independently for the COI and 16S alignments 
with MEGA 5.2 under the TN93 model, using a discrete Gamma 
(G) distribution to model differences in evolutionary rates among 
sites. The null hypotheses of equal evolutionary rates throughout 
the trees were rejected at a 5% significance level. Substitution rates 
(COI: 0.02333, 16S: 0.02528; substitutions/ma) were assigned in 
the prior settings of BEAUti 1.7.5 (part of the BEAST package), 
and uncorrelated relaxed lognormal molecular clocks were 
implemented for both sequence partitions. "Speciation: Yule 
Process" was chosen as tree prior. The BEAST analysis was 
performed with four independent runs with lO' generations each 
(sample freq.: 1000). Tracer 1.5 [34] was used to assess whether 
the four runs had converged. The four independent rvms were 
then combined with LogCombiner 1.7.5 (part of the BEAST 
package). Subsequently, 25% of the trees were discarded as burnin 
and a 50% majority rule consensus tree was calculated from the 
remaining 3x10' trees. Median node heights and 95% highest 
posterior density (95% HPD) intervals are provided for major 
nodes in the results section. The rate-calibrated linearized tree was 
then prepared for Lagrange configurator 20130526, together with 
a range matrix in which each lineage was assigned either to the 
Northern Calcareous Alps, the Southern Calcareous Alps, the 
Western Alps or the Western Carpathians. Migration was 
permitted between all regions, but lower probabilities ('0.5' instead 
of '1.0') were assigned in the dispersal constraints for migration 
between areas not being immediately adjacent (Southern Calcar- 
eous Alps and Western Carpathians; Southern Calcareous Alps 
and Western Alps; Western Carpathians and Western Alps). The 
ancestors were allowed to occupy a maximum of two geographic 



PLOS ONE I www.plosone.org 



12 



May 2014 | Volume 9 | Issue 5 | e96012 



Phylogeography of the Land Snail Orcuta dolium 




GLA1_5933 
GLA1_5934 



0.79 



SMU1_6140 
f WIW2_385_386_387_392_393_394 
A— WIW3_1161 

^^WIW3_1162 

I- WIW8_628 
r VEF1_3926 

VEF1_3929 

W VEF1_3927 
VEF 



^ VEF1 3928 




rGAI2_1108 
GAI2_1109_GAI3_1151 

1_3461 
i 1944 
1945 
3465 

— .UL4_1363_1365 
GAI3_1150_GAI4_640_642 
1364 

1509 1510 1511 
GUT3_3066_3b67 
— GUT3_3068 
GUT8_1134 
GUT8_298 
GUT8_299 
SNH21_1362 
GUT1_3041 

05 3906 
UT1_3042_3043 

G UT2_3074_3075_3076 



4A 

4B 
5 




MLF2_3919_3922_MLF3 
MLF3_3917 _3912_3913 

MLF2_3920 

MLF2_3925 
M LF2_3921_3923_MLF3_391 4 
— -H ^ MLF3_3915 
MLF6_3952 
MLF6 3953 
MLF6_3954 - 



I 0.99r "'-F2_ 

rH i~ ""-I 

0.991 "4 MLF2 
— — H I— MLi 
I J ML 
' •p mi 

^ II 



I - ML 



MLF4_3930 
MLF4_3934 

MLF4_3931 — 
MLF4_3935 
MLF4_3932 



65 



. MLF1_3938 
. MLF1_3936_3940 



1— MLF6_3944 _3941_3943 

- MLF6_3945 

— MLF6_3946 

- MLF6 3948 
MLF6_3947_3949 

- MLF6_3950 
MLF6 3951 



1A 
1B 
1C 

1D 

2 




Orcula conica 



r BGD1_4089_4090 

■ BGD2_4086 

rGAI1_1132 
GAI1_1131_1133 
r KAI1_1140 
«- KAI2_1119 

- KAI3 838 

- KWLT_1983_615_617 

■ R0F1 621 622 
< BGD2_4087 

I KAI1_1138_1982_KAI2_1120_KAI3_1107 6506 
' KWL1_614 _KA14_1142_1143J981 
• LEC1_5932 
L ROF1_620 

■ BGD7 4085 

■ BGD11_1173_BGD7_4083_4084 

■ HSC1_1392 

- MRZ6_6149 
'■MRZ5_6146 

■ HSC2_3081_MRZ6_6150_MRZ5_6147 
> YBB9 1152 
r BGD8_4126_4127_4128 
—I- BG 09.41 22 
^ BGD9_4123 

HENN10_5623 5625 ENN12 5632 5633 
ENN10_5624 _5634_ENN15_910 
ENN8_916 
' ENN8_914 
I ENN2_1135_ENN4_874_875_876_ENN6 
^ FOB3_371 _5654_ENN8_913_ENN9_5618 
r F0B4 1121 

■ GUT5_5472 

■ GUTS 5474 
- OSR1_5613 

■ SNH16 1123 

'^'^ B2_629_630_B37_FOB4_606 
5_5475 
■ OSR1_5611 
• SNH23 2814 397 
H FOB4_ll22 

rFOB8_1984 
GUT7_2936 

■ GUT6_2945 
^ GUT7_2934_2935 
h OSR3_5606 

■ FOB3_370_372_FOB5_339_FOB6_361_367 
I GUT4_5471_GUT5 5473 GUT6 2942 2943 
L OSR3_5607_5608 

fFOB8_602 
FOB8_603 
SNH23 2815 
FOB8_604 
SNH23 395 
^ SNH20_3859 
^ ENN11_893_895 
~ ENN13_1128_1129_ENN14_842 

■ ENN14_841 

- ENN15_897 

■ ENN16_1166_1167_1168 
■- ENN17_3889 

~ ENN17_3887 3890 

— ENN18_1158 
ENN18_1159 

' ENN18_1160 1 A 

- ENN19_3106 3107 

- ENN19_3108 
^ ENN1_3083_3084 

4 ENN1_3085 

ENN2_1136 
* ENN3_5646_5647_5648 

■ ENN2_1137 ENN11 894_ENN15 911 

■ ENN8_918 

■ FOB1_858 859 860 

P LAV1_865 

LAV1_864_866 
J— MRZ1_2978 
^ MRZ2 2972 2973 

- MRZ2_2974 

- MRZ4 2961 

- MRZ4_2962 

- MRZ4_2963 

- MRZ5_2970 

■ MRZ5 2971 

^ SNH11_3328_3329 
r+ NH11_3330_SNH12 3324 3325 3327 
^ SNH13_3334 _SNH13l3332l3333 

— SNH14_2947 
«— SNH15_2950 

w YBB10_3883 
HSG1_1358 

■ YBB8 95 

— SNm4_2948 

■ SNH17_315 

- SNH17_316 

■ SNH17_317 

- SNH18_1111 
r SNH19_349 

SNH19_350 
~ SNH19_351 
« HSG2_3078 

— HSC2_3079 
< SNH20_3860 

■ MRZ1_2976_2977_MRZ3_3028 

- SNH20_3862 

LSNH22_1387 
SNH22_1389 

■ SNH23_396 
^ SZK2_1520 

I BGD1_4088_BGD10_4079_4080_4081 
4 DAC1_1357_1377_1378_DAC2_1125 
I SZK2_1518_1519_SZK3 1508 SZK4 1502 1503 
I TEN3_1495 

- SZK1_1506_1507_SZK4_1501 

■ YBB10_3884 

— YBB11_1933_YBB12_1934 
P YBB12_1935 
i— YBB12 1936 

LYBB1_1148 
' YBB1_1149 
I— YBB2_1163 
^ YBB2_1165 
r- YBB4_834 
YBB5_306 308 
Tp YBB2_1 1 64_YBB3_1 126_1 127_YBB4_31 3 
YBB6_93_94 _835 YBB5 96 98 

1 YBB5_97 

■ MRZ3_3030_MRZ4_2964 

■ YBB6_92 

- YBB7 300 

■ YBB8_1124 



PLOS ONE I www.plosone.org 



13 



May 2014 | Volume 9 | Issue 5 | e96012 



Phylogeography of the Land Snail Orcula dolium 



Figure 2. Bl analysis of COi hapiotypes. The central tree provides an overview. Posterior probabilities are given at the main nodes. Black dots 
indicate nodes with maximum support. A specimen of 0. conica is used as outgroup. More details are provided in the partial trees with different 
colors corresponding to the geographic mountain regions depicted in Fig. 1. The scale bars indicate the expected number of substitutions per site 
according to the model of sequence evolution applied. 
doi:1 0.1 371/journal.pone.009601 2.g002 



areas. Alternative analyses were run with unlimited range sizes 
(aUowing the taxa to inhabit all four geographic areas) and 
assigning the same probabilities ('1.0') for migration between the 
four areas in the dispersal constraints. 

Literature Search for Fossil Records of O. dolium 

We screened various papers for fossU records of 0. dolium. The 
dating of the Alpine sites seems to be rather tentative because all 
but one site were dated only using reconciliation with vertebrate 
fossils of the same layers. Most reliable is the stratigraphic dating of 
recently investigated loess profiles of the Pannonian Basin, 
including measurements of carbonate content variations, low-field 
magnetic susceptibility and radiocarbon dating using macro 
charcoal fragments. Most of the Western Carpathian records are 
based on investigations of soil profiles as well. The publications 
with positive records of 0. dolium are listed according to the 
countries investigated: Austria [42,43], Croatia [44,45], Czech 
Republic [46-48], Germany [49], Hungary [19,50], Serbia [51- 
53] and Slovakia [47,48]. Since the evaluation of the fossil record 
of Orcula dolium was based on literature data only, no permits were 
required for that part of the study. We assigned each record to a 
time period of either cold or warm climate phases of Middle and 
Late Pleistocene according to the Quaternary divisions of the 
North European climate cycles of Zagwijn [54] : the Weichselian 
(115-11 ka ago), Saalian (350-130 ka ago) and Elsterian (475- 
370 ka ago) were considered as cold climate stages (glacials). 



compared to the warm stages (interglacials), Holocene (1 1 ka ago - 
present), Weichselian-Saalian iiiterglacial ( = Eemian; 130-115 ka 
ago) and Saalian-Elsterian interglacial ( = Holsteinian; 370-350 ka 
ago). However, it has to be noted that massive glaciations of the 
Alps and northern Europe occurred only at some times of the 
glacials, most recent during the LGM (30-18 kya). The maps 
showing the distribution of genetic clades and fossil records were 
prepared using ArcMap Desktop 10.0 and manually edited in 
Adobe Photoshop CS4 version 12. 

Results 

Mitochondrial Clades 

Among the 373 individuals of 0. dolium investigated, 197 COI 
hapiotypes are observed. The phylogenetic trees calculated with 
different algorithms show similar topologies. Ten major clades, 
some of them divided into distinct sub-clades, are highly supported 
in all analyses. Six of the major clades occur exclusively in the 
Western Carpathians (2, 3, 5, 7, 9 and 10), three clades are 
distributed solely in the Alpine region (4, 6 and 8). The final clade, 
clade 1 , is subdivided into four sub-clades themselves restricted to 
particular regions (Alps: lA; Western Carpathians: IB, IC and 
ID) (Figs. 2, 3 and 4). Sub-clade lA is distributed throughout the 
Northern Calcareous Alps but occurs also in a distinct area of the 
Southern Calcareous Alps. Similarly, clade 6 is distributed in 
separate areas of the Northern Calcareous Alps and the Southern 




Figure 3. Distribution of localities investigated in the Alps and the Western Carpathians. The symbols indicate the presence of the 
respective mt clades. If specimens of different mt clades are present at the same localities, the respective clade symbols are shown encircled with a 
line pointing towards the locality. A picture of the shell of a type specimen of 0. dolium (Syntype NHMW 14765/1820.26.61/2; size: 6.8 mm in height) 
is provided in the left upper corner. 
doi:1 0.1 371/journal.pone.009601 2.g003 
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Figure 4. Bl tree of the concatenated COI and 16S sequences. Posterior probabilities and ML bootstrap values are provided for all nodes 
above sub-clade level. The scale bar indicates the expected number of substitutions per site according to the model of sequence evolution applied. 
The dots indicate nodes with posterior probabilities of 1.0 and bootstrap values of more than 98% (black) or 95% to 97% (grey) in the ML analyses, 
respectively. 

doi:1 0.1 371 /journal.pone.009601 2.g004 



Calcareous Alps. Clade 4 also has a disjunct distribution, but in an 
east-west orientation, occurring in the eastern-most part of the 
Northern Calcareous Alps (4B) and in the Western Alps (4A). In 
contrast, clade 8 is geographically restricted to a small area in the 
Northern Calcareous Alps. The BI tree based on COI and 16S 
reveals the same highly supported clades, sub-clades and overall 
topology as the COI tree, but some nodes are better supported. 
Nevertheless, the topology remains partly ambiguous, e.g., the 
relationships between clades 1, 2 and 3 as well as between 4, 5, and 
6 are unresolved trichotomies. 

Nuclear Data Set 

The H4/H3 data is shown as a phylogenetic network (Fig. 5). It 
is not possible to incorporate the division of the mt clades into the 



H4/H3 network because the pattern differs considerably from that 
obtained in the mt trees - specimens of the same mt clades often 
provide very different H4/H3 alleles. The general scheme of the 
network exhibits three frequent alleles (HTl, HT2 and HT3), each 
encircled by several similar haplotypes, which differ by a few 
substitutions or indels only. Additionally, there are several unique 
alleles, particularly in specimens from the Western Carpathians 
(e.g. in 1996_STR2, 1376_MLF5, 3952_MLF6, 3915_MLF3 and 
1373_MLF7). In the Alpine populations, HTl is the most frequent 
allele and, apart from a unique Western Alpine allele (differing by 
2 substitutions from HTl), it is the only one found in the Western 
Alps and the Southern Calcareous Alps. In the Western 
Carpathians, HTl is detected in a single specimen (STR1_3909) 
only. However, HTl is encircled by several slightly differing 
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haplotypes of both the Northern Calcareous Alps and the Western 
Carpathians. In contrast, HT2 and HT3 are each found 
exclusively in either the Northern Calcareous Alps or the Western 
Carpathians, respectively. Thus, aU populations having HT2 
additionally feature HTl or haplotypes slightly deviating from 
HT 1 . Within the Alps, the populations of the eastern-most margin 
of the Northern Calcareous Alps exhibit the largest allele diversity. 

Genetic Diversity 

Intraspecific distances are extraordinarily high for the COI 
gene, with up to 18.4% mean /(-distance between the clades 
(Table 3). The maximum /)-distance within the Alpine populations 
is 16.9%, compared to 18.3% in the Western Carpathian ones. 
The mean /j-distance between 0. dolium and 0. conica is only slightly 
higher at 18.4% (Table 3). The respective haplotype and 
nucleotide diversities of the COI clades are high in all populations 
(Table 4). Distances were also calculated for the 16S sequences 
(Table 5). The largest mean /(-distance between clades of 0. dolium 
is 14.4%, the distance between 0. dolium and 0. conica is 19.3%, 
indicating that the conserved parts of the 1 6S suffered less from 
saturation than COI. Regarding the nc H4/H3 sequences, the 
Western Carpathian populations almost consistently show larger /i- 
distances (max.: H3 1.5%, H4 2.7%, Spacer 2.6%; clones 
excluded: H3 0.3%, H4 1.6%, Spacer 1.0%) than the Alpine 
ones (max.: H3 0.6%, H4 1.9%, Spacer 2.1%; clones excluded: 
H3 0.3%, H4 1.2%, Spacer 1.5%). The max. /(-distances within O. 
dolium (H3 1.5%, H4 2.7%, Spacer 2.6%; clones excluded: H3 
0.3%, H4 1.9%, Spacer 1.8%) are comparable to the mean 
distances between 0. dolium and O. conica (H3 0.67%, H4 1.30%, 
Spacer 2.38%). 



Molecular Clock Analysis and Reconstruction of the 
Phylogeographic Range Evolution 

The linearized tree resulting from the BEAST 1.7.5 analysis was 
combined with a biogeographical range reconstruction using 
Lagrange v.201 30526 (Fig. 6). For the nodes marking major splits 
of mt lineages, the node ages and the 95 % HPD intervals (in mya) 
are provided. The alternative ancestral subdivision/inheritance 
scenarios with likelihoods of 1 5 % or more are also indicated in the 
tree. As a main result, the analyses suggest that 0. dolium originated 
6.92 to 4.13 mya (95% HPD interval) around the Miocene- 
Pliocene boundary. However, the analyses support that the broad 
diversification into numerous lineages happened during the 
Pleistocene. According to the reconstruction of the species' 
geographic range history (ancestors allowed to occupy a maximum 
of two geographic areas/dispersal constraint for areas not 
immediately adjacent: '0.5'), the ancestral 0. dolium was distributed 
in the Western Carpathians (Maximum Likelihood of ancestral 
stage at cladogenesis event: 0.83). The alternative analysis run 
without range restrictions (ancestors allowed to occupy all four 
geographic areas/dispersal constraint for areas not immecUately 
adjacent: '0.5') predicted similar ancestral ranges but the most 
likely ancestral range scenarios generally obtained lower likelihood 
values. Moreover, alternative ancestral ranges predicted for several 
nodes comprised differing geographic areas (Fig. S2). Analyses run 
with the same range constraints as in the two previous analyses, 
but with differing dispersal constraint (migrations between all areas 
are equally likely) resulted in identical ancestral ranges for almost 
all nodes, only the ML values differed slightiy (by a maximum of 
10 percent) (data not shown). 




Figure 5. Median Joining network of the H4/H3 sequences. Bars indicate substitutions within the H3, H4 and the spacer region. The loss or 
gain of bases (indels) in the spacer region is displayed by arrows pointing in the respective direction and the numbers indicate in how many bases 
the haplotypes/alleles differ from each other. The three most common alleles are named as HTl, HT2 and HT3. The size of the circles corresponds to 
the number of sequences providing the same allele. The colors correspond to one of the mountain regions defined in Fig. 1. To facilitate the 
comparison of nc and mt data, the specimen labels and the clade affiliation are indicated next to the haplotype circles. 
doi:1 0.1 371/journal.pone.009601 2.g005 
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Table 3. Mean and maximum genetic p-distances (in %) for the COI sequences. 
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Middle and Late Pleistocene Distribution Derived from 
Fossil Record 

The information about the fossil distribution of 0. dolium is 
displayed in Fig. 7. Separate maps are presented for cold and 
warm Pleistocene climate stages because several localities provide 
records of both glacial and interglacial periods. The literature 
reports 0. dolium from about 100 fossil sites, most of which are 
located in the Western Carpathians, including 40 sites in Slovakia 
alone [47,48] . Four Czech localities provide Middle Pleistocene to 
Holocene fossils, among them the northern-most confirmed site, 
which is located about 30 km north of Prague [47] . Around 20 
sites are located in the Western Carpathians of Northern Hungary, 
with Holocene, Weichselian, Eemian and Elsterian deposits 
[19,48,50]. Some of these represent the earliest fossil records of 
the species, particularly the sites in the Hungarian Biikk mountains 
and the surroundings of Budapest. Deposits from the Slovenska 
Skala in South-Eastern Slovakia also date back to the Elsterian 
and Holsteinian about 400 ka ago [48]. Fewer data are available 
from Austria and Germany. Most of the 1 8 sites of the Northern 
Calcareous Alps are dated to the Holocene or the Weichselian. 
Only a single Alpine site (Vienna, Austria) was tentatively classified 
to the end of the Middle Pleistocene by [42]; it is not included in 
the maps because an assignment to a cold or warm climate stage is 
not possible. The Nufiloch site (Baden-Wtirttemberg, Germany) is 
the earliest fossil occurrence in the foreland of the Western Alps in 
the Early and Middle W eichselian, with an increasing frequency in 
the Late Weichselian layers [49]. Numerous fossil records of 0. 
dolium were reported from loess deposits of the Pannonian Basin 
[19,52,55], most of which are located along the Danube in 
southern Hungary and Northern Serbia, where the rivers Sava 
and Danube join. Most Pannonian sites provide Weichselian 
deposits, although the species was already present close to Novi 
Sad in Serbia at the end of the Saalian [51]. In Eastern Croatia, 
fossils of 0. dolium date back from the Weichselian to the Saalian 
(over 200 ka ago) [44,45]. 



Discussion 

Origin and Diversification of O. dolium 

Zimmermann [56] hyiJOthesized that 0. dolium originated in the 
Northern Calcareous Alps like several congeneric species. Frank 
[42] took up the same position and stated that 0. dolium emerged 
in the Northern Calcareous Alps in the Early Pleistocene. 
However, these assumptions regarding the species' origin are 
rather tentative because no Pliocene and Early Pleistocene fossils 
are known from Central Europe. The earliest record of 0. dolium 
from South-Eastern Slovakia and North- Western Hungary dates 
back to the Elsterian (475-370 kya) and even earher [19,48], 
whereas almost all records from Alpine sites were assigned to the 
WeichseUan (115-11 ka ago) or the Holocene; only a single record 
(Austria, Vienna) was vaguely dated to the Late Middle Pleistocene 
[42]. In general, the gastropod record before the Middle 
Pleistocene is extremely scarce because sediments predominandy 
consist of red clay, being inappropriate for the fossUization of 
gastropod shells [19]. Besides, environmental dynamics in the 
mountainous regions, especially in the preferred limestone 
habitats, offer few opportunities for shell fossUization. Hence, in 
the present study, hypotheses regarding the origin and diversifi- 
cation of 0. dolium are mainly based on molecular genetic data. 

Contradicting the assumption of Zimmermann [56] and Frank 
[42], the variability of the mt and nc markers and the geographic 
distribution of haplotypes support an origin of 0. dolium in the 
Western Carpathians. Seven out of ten clades occur in the 
Western Carpathians, including two highly diverged clades (9 and 
10) which split from the basal nodes of the trees. The populations 
of the Alps and the Western Carpathians are not reciprocally 
monophyletic - lineages of both areas derive from three nodes in 
the mt trees each. The number of H4/H3 alleles highly diverged 
from the three main alleles HTl, HT2 and HT3 is larger in the 
Western Carpathians as well. The geographic range reconstruc- 
tion supports a scenario in which the most recent common 
ancestor was distributed in the Western Carpathians (ML: 0.83) 
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around 6.26 mya (95% HPD: 5.15 to 7.42); the Northern 
Calcareous Alps were probably settled later (Fig. 6). A scenario in 
which the MRCA was distributed in both the Western 
Carpathians and the Northern Calcareous Alps obtained only 
low support (ML: 0.12). The distributional patterns of mt clades/ 
variants can best be explained by multiple (probably two or three) 
migrations between Alps and Western Carpathians with single 
specimens or populations carrying unique or similar mt variants, 
respectively. Since the Alpine clades are embedded within the 
diversity of the Western Carpathians, a predominant east-west 
migration route is most probable. Alternative to scenarios with 
multiple migrations, the Alpine diversity could have resulted from 
a single migration involving multiple individuals carrying strongly 
differing mt variants. However, the geographically distinct 
distributions of the Alpine mt clades suggest tiiat the lineages 
evolved independently from each other. Moreover, the molecular 
clock analysis indicated that the Alpine mt clades separated from 
their closest related Western Carpathian lineages during different 
time periods. The results of the analysis suggested that the Alpine 
mt sub-clade lA descended from the Western Carpathians rather 
recently during the Middle Pleistocene, whereas clade 8 and the 
cluster including clades 4 and 6 separated from their closest related 
Western Carpathians lineages probably during the Early Pleisto- 
cene already. The nc H4/H3 sequence patterns also support at 
least two independent migration events, as two highly diverged, 
geographically more or less separated clusters (HTl, HT2 and 
similar variants) were found. 

Pleistocene Refuges and Postglacial Expansion Routes 

One of the present study's main objectives is the detection of 
potential glacial refuges of 0. dolium. The four major limestone 
areas currently inhabited by 0. dolium (Western Alps, Northern 
Calcareous Alps, Southern Calcareous Alps and Western Car- 
pathians) are treated separately, as is the Pannonian Basin, in 
which the species is apparently not found nowadays. 

The extensive fossil record of the Western Carpathians, with 
data from both glacials and interglacials, confirms that the area 
was permanently settled, at least during the second half of the 
Middle Pleistocene and the Late Pleistocene [48] . In particular the 
extensive record of Weichselian (115-11 ka ago) fossils provides 
evidence that 0. dolium was widely distributed during the last 
glacial (Fig. 7). Moreover, despite the comparably small sample 
size, the Western Carpathian populations show a large genetic 
diversity with complex distribution patterns. There is no clear 
geographic structure regarding the distribution of mt clades and nc 
alleles, and the data do not indicate a serious loss of genetic 
diversity due to genetic bottlenecks. Unlike in the Alpine region, 
the loss of potential habitat presumably was less significant in the 
Western Carpathians, which were not affected by massive 
glaciations during Pleistocene cold stages. The scattered distribu- 
tion of limestone bedrock in the Western Carpathians is another 
factor, which may have triggered diversification and preservation 
of numerous genetic lineages. 

The eastern part of the Northern Calcareous Alps potentially 
provided the largest Alpine refuge for calciphilous taxa because it 
continuously offered non-glaciated limestone areas. Patterns of 
endemism and comparative phylogeographic analyses in Alpine 
plants [8] provide additional evidence for refuges in this area. 
Similarly, the Northern Calcareous Alps harbor a number of 
endemic species such as Trochulus oreinos and Cjlindrus obtusus which 
probably originated in that region [14]. Haplotypes of all four 
Alpine mt clades are found here, with populations located 
somewhat separated from each other, and therefore suggesting 
several smaller refuges. Moreover, the respective populations show 
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Figure 6. Reconstruction of the geographic range evolution. The map shows the distribution areas of O. dolium in the four Alpine and 
Carpathian mountain areas sampled (encoded by different colors). Small black dots represent localities sampled in the present study. The linearized 
molecular clock dated Bl tree shows the relationships of selected mt lineages (C0I/16S data) of O. dolium. Black and grey dots indicate nodes with 
high posterior probabilities (see figure for values). The colored symbols at the branch tips indicate the geographic origin of the haplotypes. The 
ancestors were allowed to occupy a maximum of two geographic areas. At the cladogenesis events (nodes), all alternative ancestral subdivision/ 
inheritance scenarios with likelihoods of 15% or more are indicated, together with the respective likelihoods, and separated by an "or". When 
scenarios for cladogenesis events involve two ancestral areas, the symbol for the likely ancestral area/— s is/are provided left to each of the two 
branches. For nodes representing major splits, node ages and 95% posterior HPD intervals are indicated. A time scale in mya is given below. 
doi:1 0.1 371 /journal.pone.009601 2.g006 
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Figure 7. Distribution of fossil O. dolium in Central Europe during warm (A) and cold (B) Pleistocene climate stages. The LGM glacier 
line indicates the nnaximum extent of glaciers during the Weichselian (35 to 1 9 ka). The grey-shaded areas correspond to the current distribution of 0. 
dolium. 
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a high diversity in nc H4/H3 alleles. The most common mt clade 
(lA) is distributed from Lower Austria (Gutenstein Alps) in the east 
to Tirol (Lechtal Alps) in the west, spanning a distance of 400 km 
(Fig. 3). A distinct mt clade (4B) is present at the eastern edge of the 
Alps (VVienerwald), a region which is geologically somewhat 
isolated from the Northern Calcareous Alps due to the predom- 
inance of sUiciclastics in the intermediate region. The large nc 
diversity with several H4/H3 alleles, each strongly diverged from 
the most common alleles HTl and HT2, and the presence of 
Weichselian fossil deposits support the assumption that the 
Wienerwald served as a long-term refug(-. The Gutenstein Alps 
at the northeastern margin of the Northern Calcareous Alps may 
have provided a refuge area for the population exhibiting mt clade 
6. Another potential refuge was located close to the former LGM 
glacier line in Salzburg and Upper Austria - mt clade 8 is 
distributed exclusively in this region, and fossil deposits indicate 
the local presence of the species during the Weichselian [42]. 
Some of the specimens carrying haplotypes of mt clade 8 possess 
the nc allele HT2 or similar ones, which elsewhere occur in the 
Wienerwald only. This can either be explained by past gene flow 
between the two currently separated populations or by ancestral 
polymorphism, i.e., the persistence of ancestral histone variants in 
both areas. 

The geographic range reconstruction suggests that the Western 
Alps were settled from the Northern Calcareous Alps during the 
Middle Pleistocene 1.05 to 0.66 mya (95% HPD interval). Samples 
from the Western Alpine sites all form a single mt sub-clade (4A) 
and have the nc haplotype HTl or similar types. Although 
material is available from three sample sites only, the Western 
Alpine populations show higher distances in the mt sequences 
(max. uncorrected jft-distance COl: 3.8%) than any other Alpine 
clade. The unique presence of the highly diverged mt clade 4A 
provides support for a refuge in the Western Alps. Since the 
Western Alps were almost completely covered by ice during 
several glacials, populations may have outlasted the glacial periods 
in several smaller refuges at the Western Alpine margins of 
Switzerland and France as was proposed for Trochulus villo.ms [15] 
or Carychium tridentatum [13]. Fossils in Early to Late Weichselian 
deposits of the Western Alpine foreland (NuPloch, Baden- 
Wiirttemberg, Germany) clearly support that assumption, at least 
for the last glacial period [49]. We have no molecular data from 
the German and French areas, but a common ancestry of the 
western populations is supported by similar conchological traits 
(collection material of the Natural History Museum, Vienna and 
Naturmuseum Senckenberg, Frankfurt am Main; Harl et al. in 
prep.). 

The Southern Calcareous Alps were almost completely covered 
by glaciers during the LGM and hence provided only a small 
potential refuge for calciphilous taxa in the eastern-most part: a 
small non-glaciated area in the Karawanks [57]. The phylogenetic 
data argue against a permanent refuge of 0. dolium in the Southern 
Calcareous Alps and suggest that clades lA and 6 probably 
reached this region rather recently. In the Southern Calcareous 
Alps, sub-clade lA is found at a single site (GAIl) only, with 
haplotypes similar to those in populations of the Tirolian and 
Bavarian Northern Calcareous Alps. Since both areas were fully 
glaciated, clade lA might have crossed the Central Alps in this 
area after the LGM. Apart from this single locality, all other sites 
of the Southern Calcareous Alps possessed haplotypes of clade 6, 
which is distributed in the Northern Calcareous Alps as well. The 
similarity of clade 6 haplotypes rather indicates a very recent 
expansion during the Late Pleistocene or Holocene. However, the 
presence of a distinct variant of clade 6 at a single site in the 



Karawanks (KWNl) could be an indication for a Southern 
Calcareous Alpine refuge. 

The fossil record indicates a more or less continuous presence of 
0. dolium in the Pannonian Basin at least from the Saalian onwards 
(over 200 ka ago). The habitats of the Pannonian populations were 
probably patchily distributed forests near rivers, which were 
present in the area even during the LGM [7] . The occurrence of 
trees at the respective sites is additionally supported by the co- 
occurrence of other woodland species in the same loess strata, such 
as Semilimax semilimax, Ena montana or Aegopinella ressmanni [45,58]. 
One might ask \vhether riparian drift from Alpine or Western 
Carpathian regions could account for the presence of 0. dolium in 
the fossil record of the southern part of the Pannonian Basin. 
However, the high abundance of fossil 0. dolium in the Pannonian 
Basin indicates a local source. The contemporary absence of the 
species is probably the result of anthropogenically induced loss of 
suitable habitat. Lozek [48] stated that deforestation and 
dehydration are probably the reasons why thes(; areas lack several 
gastropod species which were still widely distributed during the 
Eemian. Thus, the expansion of agricultural areas is a reasonable 
explanation for the decline of O. dolium populations in the 
Pannonian Basin during the Holocene. 

Genetic Differentiation and Taxonomic Considerations 

The intraspecilic distances measured for the mt genes are 
among the highest found in pulmonate species (uncorrected p- 
distances: COI, 18.4%; 16S, 14.4%). By comparison, the 
genetically diverse hehcid taxa Theba pisana and Arianta arbustomm 
show COI distances of 13.6% [59] and 15% [10], respectively. 
Regarding the non-protein coding 16S, higher intraspecilic 
distances were found in the clausiliid species Charpentieria itala with 
24.2% [60]. 16S divergence is also high in the bradybaenid 
Euhadra quaedta with 14.1% [61]. Regarding the nc sequences 
analyzed, the largest ^-distance measured within the protein 
coding H4 and H3 sequences is 0.8%. The highest distance 
observed in the non-protein coding spacer region is 1.8%. [59] 
reported jS-distances of 0.5% in the non-coding ITSl sequences of 
Thebn pisniia, \vhereas the intraspecific sequence divergence within 
Arion subfusfus, a species extremely variable in its mtDNA, is only 
0.3%. for the ITSl sequence [62]. 

Considering the large genetic variability found within the 
populations of 0. dolium, the question arises whether some of the 
hneages might even represent distinct species. However, there are 
no indications of hybridization barriers that would suggest sphtting 
the groups into different species. The data indicate g(;ne flow 
between clades, as suggested by the fact that the specimens 
displaying the main histone gene variants each feature haplotypes 
of very distant mt clades (see Fig. 5). However, whether genetic 
groups correspond to currently accepted subspecies remains a 
problematic issue. More than 20 subspecies have been described 
for 0. dolium, equally divided between the Alps and the Western 
Carpathians [20]. Most were characterized by minor differences in 
shell shape and the formation of the aperture folds. Our study 
includes specimens from several type localities of Slovakian 
subspecies, namely of 0. d. titan Brancsik, 1887, 0. d. brancsikii 
Clessin, 1887, 0. d. minima Brancsik, 1887, and 0. d. cebratica 
Westerlund, 1887 (Table 1 and Fig. SI). However, none of the 
clades can be definitively attributed to one of these subspecies. For 
instance, populations of the very slender, large-shelled 0. d. 
brancsikii share the same mt haplotypes (clade IC), and the nc 
haplotype HT3, with compact, small 0. d. minima morphs. The 
very large-shelled 0. d. titan possesses a diverged mt haplotype 
within clade 9 but shows the most common nc haplotype HTl, 
which otherwise is found in the Alps only. The Slovakian 0. d. 
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cebratica (clade 5) clusters with the Alpine mt clades 4 and 6 but 
displays a distinct nc haplotype (HT3). Despite the extremely high 
conchological and genetic variability in the Western Carpathians, 
we could not detect a clear correspondence between genetic 
haplotype groups and subspecies ranges defined by [17]. The 
geographically isolated and conchologically aberrant endemic of 
the Wienerwald, 0. dolium infirm, is genetically differentiated fi-om 
other populations in the mt trees (sub-clade 4B). In contrast, 
individuals from the type localities of 0. dolium gracilior Zimmer- 
mann, 1932, 0. dolium pseudogularis WsLgner, 1912, 0. dolium edita 
Ehrmann, 1933, and 0. dolium raxae Gittenberger, 1978, all possess 
haplotypes of the homogeneous mt clade 1 and exhibit the nc 
haplotypes HT 1 or its derivatives, none of them forming a distinct 
sub-clade. The latter two taxa, 0. d. edita and 0. d. raxae, initially 
were of special interest for our study because they were reported to 
occur only at high altitudes [17, ,56]. Nonetheless, specimens from 
the corresponding localities do not differ from the common types 
of the surrounding lowlands in the markers analyzed. However, 
for final taxonomic decisions, the most important question is 
whether the presumed morphological distinctness of the various 
subspecies can be confirmed by morphometric investigations. 
Such analyses are under way to quantify morphological differences 
and to focus on the effect of altitude on shape formation and shell 
size (Harl et al., in prep.). 

Supporting Information 

Figure SI Pictures of selected types of O. dolium 
subspecies. Specimens collected at the type localities of the 
respective subspecies were investigated in the present study. The 
types shown shall rather be considered as examples for the species' 
variability than defined discrete entities of morphologically 
separated populations. Some of the types indeed represent extreme 
morphs but transitional forms are found in most populations. The 
pictures were already published by Harl et al. (201 1) together with 
data on all other currently known subspecies. In the following we 
provide the collection data of the specimens show n: A: 0. d. dolium 
(syntype NHMW 14765/1820.26.61/2), B: 0. d edita (syntype 
LML ALT/5319/1), C: 0. d raxae (syntype LML ALT/5354/1), 
D: 0. d. pseudogularis (syntype NHMW 56158), E: 0. d. gracilior 
(syntype LML ALT/5343), F: 0. d. infima (syntype LML/ 
ALT5353/1), G: 0. d brancsikn (syntype ? NHIv/nV J. N. 22075), 
H: 0. d. titan (syntype NHMW 68377 (5448)/3), I: 0. d cebratica 
(syntype MNHG Wsdd2090), J: 0. d minima (syntype 27044/2). 
Abbreviations for Museums: NHMW (Naturhistorisches Museum, 
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Wien), LML (Oberosterreichisches Landesmuseum, Linz) and 
MNHG (Naturhistoriska Museum, Goteborg). The scale bar 
indicates 5 mm. 
(TIF) 

Figure S2 Reconstruction of the geographic range 
evolution. The map shows the distribution areas of 0. dolium in 
the four Alpine and Carpathian mountain areas sampled (encoded 
by different colors). Small black dots represent localities sampled in 
the present study. The linearized molecular clock dated BI tree 
shows the relationships of selected mt lineages (COI/16S data) of 
0. dolium. Black and grey dots indicate nodes with high posterior 
probabilities (see figure for \'alues). The colored symbols at the 
branch tips indicate the geographic origin of the haplotypes. The 
ancestors were allowed to occupy ail four geographic areas. At the 
cladogenesis events (nodes), all alternative ancestral subdivision/ 
inheritance scenarios with likelihoods of 15% or more are 
indicated, together with the respective likelihoods, and separated 
by an "or". When scenarios for cladogenesis events invok e two or 
more ancestral areas, the symbol for the likely ancestral ar<;a/ — s 
is/are provided left to each of the two branches. For nodes 
representing major sphts, node ages and 95% posterior HPD 
intervals are indicated. A time scale in mya is given below. 
(TIF) 
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